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£ I ABSTRACT 

C$ . Chemical evolution models are presented for the anomalous globular cluster to Centauri. After 

demonstrating that the chemical features of u> Cen can not be reproduced in the framework 
1"*** \ of the closed-box self-enrichment scenario, we discuss a model in which this cluster is the 

remnant of a dwarf spheroidal galaxy evolved in isolation and then swallowed by the Milky 
Way. Both infall of primordial matter and metal-enriched gas outflows have to be considered 
in order to reproduce the stellar metallicity distribution function, the age-metallicity relation 
and several abundance ratios. Yet, as long as an ordinary stellar mass function and standard 
stellar yields are assumed, we fail by far to get the enormous helium enhancement required 
to explain the blue main sequence (and, perhaps, the extreme horizontal branch) stellar data. 

■ Rotating models of massive stars producing stellar winds with large helium excesses at low 
\ metallicities have been put forward as promising candidates to solve the 'helium enigma' of 

oj Cen (Maeder & Meynet, 2006, A&A, 448, L37). However, we show that for any reasonable 
choice of the initial mass function the helium-to-metal enrichment of the integrated stellar 
population is unavoidably much lower than 70 and conclude that the issue of the helium 
enhancement in oj Cen still waits for a satisfactory explanation. We briefly speculate upon 
O ' possible solutions. 

r/3 , Keywords: galaxies: dwarf- galaxies: evolution -globular clusters: individual (u> Centauri) 

■ - stars: abundances - stars: chemically peculiar. 
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■ 1 INTRODUCTION sodium, magnesium and aluminum could be due to proton-capture 
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fusion reactions occurring during quiescent hydrogen and helium 

burning within the giants themselves, variations in the abundances 
dow into astrophysics (see Smith 2004, for a recent review). The r , , . . , , 

r J 6 ol heavier species - such as iron-peak and neutron-capture ele- 

ments - should be immune to such processes and, therefore, likely 
reflect patterns imprinted on the observed stars by previous stellar 
generations (Lloyd Evans 1983; Cohen & Bell 1986; Norris & Da 
Costa 1995; see also Gratton, Sneden & Carretta 2004, for a com- 
prehensive review on abundance variations within GCs). 



The globular cluster (GC) wCen (NGC5139) is a unique win- 

estimated mass, between (2.5 ± 0.3) x 10 6 M Q (van de Ven 
et al. 2006) and 5.1 x 10 6 M (Meylan et al. 1995) - or even 
7.1 x 10 6 M Q (Richer et al. 1991), makes it the most massive 
GC of the Milky Way. Its total mass better compares with that 
of a small dwarf spheroidal galaxy such as Sculptor (^sculptor — 
6.4 x 10 6 M Q ; Mateo 1998) rather than with that of a typical 

Galactic GC (^#ggc ~ 10 s M ; Harris 1996). The large degree The first indication that u Cen was chemically inhomogeneous 

of chemical self-enrichment observed in u Cen giant and subgiant cam e from the large color width of the red giant branch (RGB) es- 

members also sets it apart from all the other Galactic globulars. tablished in an earlier photometric work by Woolley (1966). In the 

Indeed, w Cen is the only known GGC to exhibit a large degree last decade, both extensive spectroscopic surveys and wide-field 

of chemical self-enrichment in all the elements studied (Freeman & photometric studies - mostly of giant u Cen members - have defi- 

Rodgers 1975; Cohen 1981; Mallia & Pagel 1981; Gratton 1982; nitely confirmed the existence of several (up to five) discrete stellar 

Norris & Da Costa 1995; Smith, Cunha & Lambert 1995; Smith populations covering a large range in metallicity (Norris, Freeman 

et al. 2000; Pancino et al. 2002). While star-to-star variations in & Mighell 1996; Suntzeff & Kraft 1996; Lee et al. 1999; Hilker 

the abundances of light elements such as carbon, nitrogen, oxygen, & Rkhtler 2000; Hughes & Wallerstein 2000; Pancino et al. 2000; 

Rey et al. 2004; Sollima et al. 2005a). The metallicity distributions 
of (sub)giant stars and main sequence tumoff (MSTO) stars look 

* E-mail: donatella.romano@oabo.inaf.it pretty much the same (Stanford et al. 2006). From their spectro- 
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scopic and photometric MSTO data, Stanford et al. (2006) find that 
the formation of u> Cen took place most likely over 2-4 Gyr; both 
a null age range and age ranges higher than 6 Gyr are deemed un- 
likely (cf. previous works by Norris & Da Costa 1995; Hilker & 
Richtler 2000; Hughes & Wallerstein 2000; Smith et al. 2000; Pan- 
cino et al. 2002; Rey et al. 2004; Sollima et al. 2005a; see also 
Kayser et al. 2006). It has also been found that the faint main se- 
quence of ui Cen splits into at least two distinct branches (Ander- 
son 1997; Bedin et al. 2004; Sollima et al. 2006b). The blue main 
sequence (bMS) contains ~25 per cent of the stars and is 0.3 ± 
0.2 dex more metal-rich than the red one (rMS; Piotto et al. 2005). 
These observations are most likely explained by an anomalously 
high helium abundance of bMS stars of Y > 0.38 (AY/AZ ^ 70; 
Norris 2004; Piotto et al. 2005). Such an intriguing possibility has 
been quantitatively investigated in the framework of somewhat ide- 
alized scenarios for the formation and evolution of ui Cen (Bekki & 
Norris 2006). 

It has been shown (Norris et al. 1997) that the metal-rich com- 
ponent of the cluster, which is more centrally concentrated, has a 
smaller line-of-sight velocity dispersion and a lower systemic ro- 
tation about the cluster's minor axis than the metal-poor one. A 
more recent analysis by Sollima et al. (2005b) confirms the trend 
of decreasing velocity dispersion with increasing metal abundance 
in the metallicity range —2.0 < [Fe/H] < —1.0, but shows also that 
in the extreme metal-rich extension of the stellar population, now 
better sampled thanks to new generation instrumentation (Pancino 
et al. 2000), this decreasing trend is reversed (see fig. 9b of Sol- 
lima et al. 2005b). Asymmetries in the distribution and velocity of 
the stars could testify past accretion events within to Cen (Ferraro, 
Bellazzini & Pancino 2002; Pancino et al. 2003), but evidence for 
these is not definitive yet (Platais et al. 2003). To further compli- 
cate the overall picture, the orbit of ui Cen is found to be strongly 
retrograde, almost coplanar with the Milky Way disc, and to have 
small apogalacticon, unlike any known Galactic GC (e.g. Majewski 
et al. 2000). At variance with most globulars, the relaxation time for 
a; Cen is very long, up to a few times 10 9 years in the core and a 
few times 10 10 years at half-mass radius (Meylan et al. 1995; van 
de Ven et al. 2006). Indeed, recent results suggest that the cluster 
is not yet relaxed even in the central regions (Ferraro et al. 2006). 
This cluster could then mantain for a fairly long time the imprinting 
of its initial conditions, thus allowing one to use the currently ob- 
served distribution of stars of different populations in order to trace 
back the cluster formation and evolution (e.g. Merritt, Meylan & 
Mayor 1997). 

The large mass, spread in element abundances, flattened shape 
and rotation all come close to the picture where to Cen is the sur- 
viving remnant of a larger system. The chemical and kinematical 
segregations detected in uj Cen add even more relevance to the pic- 
ture of a dwarf galaxy progenitor being subject to accretion events, 
since such gradients are present in nearly all of the Local Group 
dwarf spheroidals (e.g. Harbeck et al. 2001; Tolstoy et al. 2004; 
Koch et al. 2006, to name a few). Several authors have speculated 
on the possibility that u> Cen is the naked nucleus of a dwarf satel- 
lite galaxy captured into a retrograde Galactic orbit many billion 
years ago (Dinescu, Girard & van Altena 1999; Majewski et al. 
2000; Smith et al. 2000; Gnedin et al. 2002; Bekki & Norris 2006), 
following more general ideas on an accreted origin for GCs in our 
own as well as external galaxies (Zinnecker et al. 1988; Freeman 
1993). Self-consistent dynamical models as well as N-body hydro- 
dynamical simulations can be found in the literature which succeed 
to reproduce the main features of the cluster by assuming that it 
formed in isolation and then fell inside the Galactic potential well 



(e.g. Carraro & Lia 2000; Bekki & Freeman 2003; Tsuchiya, Kor- 
chagin & Dinescu 2004). The total mass of the parent object ranges 
from 10 s Mq to some 10 9 M@ in those models. Alternative sce- 
narios envisage an off-centre stellar supercluster seed, which would 
have trapped older galactic field stars during its formation process 
(Fellhauer, Kroupa & Evans 2006). Other possible explanations for 
the cluster origin - namely, merging between two or more smaller 
globulars or between a dwarf galaxy and a cluster (Icke & Alcaino 
1988; Norris et al. 1997), or formation triggered by cloud-cloud 
collisions (Tsujimoto & Shigeyama 2003) - have been put for- 
ward, but seem less likely on the basis of the continuous trends 
of heavy elements-to-iron and lanthanum-to-iron observed for S 
stars in a; Cen (Vanture, Wallerstein & Suntzeff 2002, and refer- 
ences therein). Clear evidence against the parent system evolving 
as a closed box has been put forward by Ikuta & Arimoto (2000). 

In this paper we deal with the chemical evolution of ui Cen. We 
discuss two possible scenarios of formation: in the first one, the pre- 
cursor of uj Cen is a small system which evolves either as a closed- 
box or with some exchange of matter with the surroundings; in the 
second one, the cluster is the leftover of an ancient nucleated dwarf 
galaxy swallowed by our Galaxy some 10 Gyr ago. We pay special 
attention to the intriguing subject of helium enhancement in bMS 
stars. For the first time, this critical issue is faced in the framework 
of a complete, self-consistent chemical evolution model, taking all 
the relevant physics into account. In Section 2 we list the observa- 
tions that we use to constrain our chemical evolution model. The 
model is presented in Section 3. Section 4 summarizes the model 
results, that are discussed and compared to previous investigations 
in Section 5. In Section 5 we also draw our conclusions. 



2 THE NATURE OF THE CHEMICAL ENRICHMENT IN 
OMEGA CENTAURI 

It was recognized several years ago that simple models for clus- 
ter enrichment can not reproduce the number of metal-rich stars 
observed in the metallicity distribution function (MDF) of ui Cen 
giants (Norris et al. 1996; Suntzeff & Kraft 1996; Ikuta & Arimoto 
2000). Based on Ca abundances obtained from low-resolution spec- 
tra, Norris et al. (1996) found a bimodal distribution with a metal- 
poor component, peaking at [Ca/H] ~ — 1 .4 dex and comprising 
nearly 80 per cent of the stars, and a metal-rich one, comprising 
nearly 20 per cent of the stars at [Ca/H] ~ —0.9 dex (Fig. 1, lower 
right panel). Suntzeff & Kraft (1996) also found a sharp rise at low 
metallicities and a high-metallicity tail, but no evidence for a sec- 
ondary hump at higher metallicities (Fig. 1, upper right panel). 

The sharp rise to a mean of [Fe/H] ~ — 1.6 dex and the long 
tail at higher metallicities have been confirmed by a number of sub- 
sequent studies, using both higher resolution spectroscopy and/or 
high quality photometry (e.g. Hilker & Richtler 2000; Frinchaboy 
et al. 2002; Sollima et al. 2005a; Stanford et al. 2006; Kayser et 
al. 2006). Some of these works also show that the MDF is more 
complex than previously thought, with several separate peaks iden- 
tified in the observed distribution (e.g. Sollima et al. 2005a; Fig. 1, 
left panel). According to subgiant branch (SGB) data, the differ- 
ent populations have ages comparable within 2 Gyr; actually, they 
might be even coeval (Sollima et al. 2005b; see also Ferraro et al. 
2004). A wider age range, At ~ 2-4 Gyr, is inferred from MSTO 
data (Stanford et al. 2006; see also Hilker et al. 2004 and Kayser et 
al. 2006). 

The abundance ratios of different chemical species provide an- 
other independent constraint on the evolutionary time-scales. As 



© 2006 RAS, MNRAS 000, 1-1 1 



The chemical evolution ofuj Cen 3 



[M/H] 



[Fe/H] 

-1.5 -1 -0.5 



Sallima et al. 2005a -Suntzeff & Kraft 1996- 
(1364 red giants) _l (199 subgiants) 



n. 





Norris et al. 1996 
(517 red giants) 



0.4 

0.35 

0.3 

0.25 

0.2 

0. 1 5 

0.1 

0.05 

0.--I 
0.35 
3 
0.35 
2 
0.15 
0.1 
0.05 



-1.5 -1 



[M/H] 



-1.5 -1 -0.5 



[Ca/H] 



Figure 1. Shown are some of the observed metallicity distribution func- 
tions for lu Cen giants. General features of all the observed distributions are 
the steep rise at low metallicities and the high-metallicity tail (see text for 
details). 



first discussed by Lloyd Evans (1983), the enrichment of .s-process 
elements (e.g. Ba, La) in the metal-rich lu Cen stars is likely to 
be that of the gas clouds out of which the stars formed. Smith et 
al. (2000) suggested previous generations of low-mass asymptotic 
giant branch (AGB) stars as polluters and argued in favour of a 
protracted period of star formation in a; Cen, of the order of 2- 
3 Gyr. Shorter time-scales are inferred from the observed knee in 
the [a/Fe] versus [Fe/H] relation, if Type la supernovae (SNela) in 
lu Cen restore the bulk of their iron to the interstellar medium (ISM) 
on the same time-scale as in the solar neighbourhood (~ 1 Gyr; Pan- 
cino et al. 2002). However, one must be aware that the time-scale 
for the maximum enrichment by SNela in a specific system de- 
pends strongly on the assumptions about the SN progenitors, stel- 
lar lifetimes, initial mass function (IMF), star formation rate (SFR) 
and, last but not least, possible metal-enriched gas outflows from 
the system. It has been demonstrated (Matteucci & Recchi 2001, 
and references therein) that time-scales as long as ~4— 5 Gyr can 
be obtained with suitable assumptions. We will come back to this 
issue later on, when discussing model results in Section 5. 

Gnedin et al. (2002) used simple dynamical modeling to 
demonstrate that if lu Cen had always evolved in isolation on its 
present orbit in the Milky Way, it would have been unable to retain 
the i'-process-rich wind material from AGB stars because of strip- 
ping from ram pressure during many passages through the disc. 
Subsequent numerical simulations demonstrated that an lu Cen-like 
object can originate from a nucleated dwarf galaxy intruding into 
the Milky Way. Both the orbital parameters and the observed sur- 
face brightness profile of the present-day lu Cen are reproduced by 
a tidal disruption scenario where the falling dwarf has its outer stel- 
lar envelope almost completely stripped, whereas a central, dense 
nucleus still survives owing to its compactness (Bekki & Freeman 
2003; Tsuchiya et al. 2004; Ideta & Makino 2004). Interestingly 
enough, also the number fraction (25 per cent) of the bMS stars 
with extreme helium enhancement (AY ~ 0. 1 2-0. 14; Norris 2004; 



Piotto et al. 2005) can not be explained as long as they originated 
from the ejecta of rMS stars initially within lu Cen. Rather, most of 
the helium-rich gas necessary to form the bMS had to come from 
field stellar populations surrounding lu Cen when it was the com- 
pact nucleus of a Galactic dwarf satellite (Bekki & Norris 2006). 
In other words, the original total mass of the rMS population must 
have been larger than the present-day one and most of the rMS stars 
must have been removed from the proto-oj Cen after their ejecta 
were used to form the bMS population. Searches for tidal debris 
from uCen's hypothetical parent galaxy in the solar neighbour- 
hood are possible in principle (Dinescu 2002). Indeed, a distinct 
population of stars with lu Cen-like phase-space characteristics and 
metallicities consistent with those of lu Cen members emerges from 
catalogs of metal-deficient stars in the vicinity of the Sun (Dinescu 
2002; Meza et al. 2005). 

Though gas fueling from an ancient host galaxy stands as an 
attractive hypothesis in many respects, to the best of our knowl- 
edge this possibility has never been studied by means of fully self- 
consistent chemical evoution models before. 



3 THE CHEMICAL EVOLUTION MODEL 

We follow the evolution of the abundances of several chemical 
species in the gaseous medium out of which lu Cen's stars form, 
in the case of either a GC precursor or a dwarf spheroidal pro- 
genitor. The adopted chemical evolution model is one zone, with 
instantaneous ^ and complete mixing of gas inside it and no instan- 
taneous recycling approximation (i.e., the stellar lifetimes are taken 
into account in detail). The GC precursor has an initial mass of 
the same order of magnitude of the present one (J'uan = 2.5- 
5 x 10 6 Mq). Both a closed-box model and models where ex- 
changes of matter with the surroundings are allowed are analyzed. 
For the dwarf galaxy precursor, when the computation starts the 
baryonic matter (cooling gas) is embedded in a relatively massive 
(^dark/^bai = 10), diffuse (i?d;uk/^cff = 10), virialized dark mat- 
ter halo. Initially, the baryonic mass inside the dark matter poten- 
tial well is two orders of magnitude higher than the present one. 
As soon as the star formation begins, the thermal energy of the gas 
starts to increase as a consequence of multiple SN explosions, even- 
tually exceeding its binding energy^ . When this condition is met, 
part of the gas escapes the galactic potential well; we assume that 
this gas is definitively lost from the system. 



3.1 Basic equations 

We use the following basic equation (Tinsley 1980) 



dt 



-Xi(t)iP(t) + m(t) 



dt 



dt 



(1) 



to track the evolution of the fractional gas mass in the form 
of element i normalized to the initial gaseous mass, Sfj(i) = 
Xi(t)^#gas(t)/^bar- The quantity Xi(t) represents the abundance 



t Notice that 'instantaneous ' actually means 'on time-scales shorter than 
the adopted timestep for integration of the equations'. 
t We compute the binding and thermal energies of the gas according to the 
recipes of Bradamante, Matteucci & D'Ercole (1998). However, we adopt 
a typical efficiency of thermalization from both Type II and Type la SNe of 
??SNII = ^SNla = 0.20, rather than 0.03 (see Romano, Tosi & Matteucci 2006, 
and references therein). 
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by mass of the element i at the time t; by definition, the sum- 
mation over all the elements in the gas mixture is equal to unity. 
The first term on the right hand side accounts for gas consumption 
by star formation; the second term on the right hand side refers 
to gas return by dying stars. All the complicated dependencies on 
the adopted stellar initial mass function and lifetimes, SNIa pro- 
genitors and stellar nucleosynthesis products hidden in the Mi it) 
term are not made explicit here; the interested reader can find all 
of them discussed in considerable detail in Matteucci & Greggio 
(1986). Suffice it here to say that in this work we use an extrapo- 
lated Salpeter (1955) IMF or a Scalo (1986) IMF, both normalized 
to unity over the 0.1-100 M stellar mass range, to show the pre- 
dictions from both steep and flat, although 'standard', IMFs. The 
adopted stellar nucleosynthesis prescriptions are discussed in Sec- 
tion 3.3. 

In our models the SFR is a simple Schmidt's (1963) law: 

^(t) = v<3 k {t). (2) 

The quantity v is the star formation efficiency, namely the inverse 
of the typical time-scale for star formation, and is expressed in units 
of Gyr -1 . The exponent k is set to be 1. 

The last two terms in Equation (1) account for any gas inflow 
and/or outflow. For the closed system, all the gas available for star 
formation is in situ when the star formation begins at t = 0, and no 
infall from outside is considered. For the open systems, the rate of 
gas infall is parametrized as 



<JS^(t) 



X? 



-t/T 



dt 



T(l 



V/T) ■ 



(3) 



with r, the infall time-scale, set to be 0.5 Gyr, and X™, the abun- 
dances of the infalling gas, set to their primordial values. In par- 
ticular, we assume Y P = 0.248, in agreement with the predictions 
of the standard big bang nucleosynthesis (SBBN) theory and the 
constraints from the cosmic microwave background (CMB; see Ro- 
mano et al. 2003, and references therein). 

The rate of gas loss via SN-driven large-scale outflows is dif- 
ferent from zero only for the open models, and simply proportional 
to the amount of gas present at the time t: 



dt 



= WiXi(t)&{t). 



(4) 



The quantity Wi is a free parameter which describes the effi- 
ciency of the galactic wind; it is expressed in Gyr -1 and may have 
different values for different elements (e.g. Recchi, Matteucci & 
D'Ercole 2001). 



3.2 The evolutive context 

According to Bekki & Freeman (2003, and references therein), in 
the dwarf progenitor scenario the initial stellar mass of lo Cen's host 
was significantly higher than that currently observed and can be 
estimated as 



(5) 



(1 — /lost)/n ' 



where /„ = 0.05 is the mass fraction of the compact nucleus 
and /i os t = 0.2 is the stellar mass fraction that gets lost through 
long-term (~10 Gyr) tidal interaction with the Milky Way. If 
^Lccn — 5 x 10 6 M Q (Meylan et al. 1995), one gets ^#dwarf = 
1.25 x 10 8 Mq. Since cen might be a factor of two lower (van 
de Ven et al. 2006), we run several models, starting with an initial 
gaseous mass J(^ s (t = 0) = ./^bar = 5 x 10 s — 10 9 M Q and ending 



up with a stellar mass ^Ctars(i = 3) ~ 5 x 10 7 -10 8 M . Given 
the similarity of the results, in the case of the dwarf galaxy parent 
we will show only the predictions for our most massive model (y = 
0.35 Gyr" 1 , w™ x ^Gyr" 1 ). 

Though in our models the star formation activity lasts ~3 Gyr, 
it is worth noting that most of the stars (nearly 75 per cent of the 
cluster population) actually form during the first 1 Gyr. In the dwarf 
progenitor scenario, while the star formation proceeds, the galaxy 
gets almost completely depleted of its gas, partly owing to gas con- 
sumption by the adopted long-lasting star formation activity, but 
most of all because of efficient gas removal through the large-scale 
galactic outflows. 

The chemical properties of our ui Cen progenitor systems are 
discussed and compared with the available observations in Sec- 
tion 4. 



3.3 Nucleosynthesis prescriptions 

In this work we adopt the metallicity-dependent yields of van den 
Hoek & Groenewegen (1997) for single low- and intermediate- 
mass stars (LIMSs; 0.9 ^ 7ti/Mq 8) and the yields of Nomoto 
et al. (1997) for massive stars (13 < m/M < 70). These latter are 
computed for a solar chemical composition of the stars. The stel- 
lar yields are then scaled to the current metallicity of the model 
by means of the production matrix formalism (Talbot & Arnett 
1973). We (arbitrarily) use linear interpolations and extrapolations 
to cover the 1-100 M stellar mass range. The effect of adopt- 
ing different yield sets will be thoroughly analyzed in a forthcom- 
ing paper (Romano et al., in preparation; notice that changing the 
adopted yield sets is not expected to affect significantly the main 
results presented in this paper). 

For stars evolving in binary systems which will give rise 
to SNIa explosions, the nucleosynthesis prescriptions are from 
Iwamoto et al. (1999). In our models, a substantial fraction of Cu 
and Zn is produced by SNela, following the suggestions of Mat- 
teucci etal. (1993). 



4 RESULTS 

4.1 The closed-box picture 

In this section we briefly discuss the results obtained in the frame- 
work of the closed-box self-enrichment scenario, where any matter 
exchange between the proto-o; Cen and its environment is strictly 
forbidden. 

In Fig. 2, the theoretical metallicity distribution as a function 
of [Fe/H] or [Ca/H] (thick solid lines) is compared with the photo- 
metric and spectroscopic empirical ones (thin dashed histograms; 
Norris et al. 1996; Suntzeff & Kraft 1996; Sollima et al. 2005a). 
This comparison shows the shortcomings of this simple picture. 
In the context of a closed-box evolution, metallicities much higher 
than observed are quickly attained, so that most of the stars form 
from matter with a chemical composition from one tenth of solar to 
supersolar, at variance with the observations. 

The results displayed in Fig. 2 are for the Scalo (1986) IMF, 
that - owing to its steeper slope for m > 2 M with respect to 
Salpeter's, x = 1.7 rather than 1.35 - allows a lower fraction of 
high-mass stars to pollute the ISM with their metal-rich ejecta. Ob- 
viously, by adopting a Salpeter IMF the predicted distribution goes 
towards even higher metallicities. Flattening the IMF in the very 
low stellar mass domain, as in a Kroupa et al.'s (1993) IMF, does 
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Figure 2. Comparison of observed (thin dashed histograms) and predicted 
(thick solid lines) metallicity distribution functions of ui Cen stars. Theoret- 
ical predictions refer to the closed-box self-enrichment scenario. The the- 
oretical distributions are convolved with a Gaussian of dispersion cr = 0.2 
dex in order to (generously) take observational errors into account. 

not prevent the system from suddenly reaching too much high metal 
abundances. Even in the framework of extreme scenarios where the 
formation of SNIa progenitors is suppressed and lower Fe yields 
from core collapse SNe are assumed (Chieffi & Limongi, in prepa- 
ration), we still fail to reproduce the observed MDF. Since within 
the scheme of a closed-box self-enrichment we can not fulfil even 
such a basic observational constraint, we deem it meaningless to 
analyze further the model results and switch to the 'open' scenar- 
ios. 

Indeed, only if a substantial fraction of SN ejecta escapes from 
the cluster, the observed MDF can be reproduced. Our result sup- 
ports the findings of Ikuta & Arimoto (2000) that significant out- 
flow from the cluster is needed in order to reduce the effective yield 
per stellar generation. Ikuta & Arimoto (2000) obtained the best fit 
to the observed MDF with a model involving gas outflow, infall 
at the very early stage of chemical evolution and a bimodal IMF. 
However, the duration of star formation for their best-fit model was 
only 0.28 Gyr, much shorter than that inferred from current obser- 
vations and assumed here. 

4.2 A stripped dwarf galaxy? 

An open, small-mass model hence offers a viable solution, but the 
questionable assumption has to be made that, while SN products 
easily leave the cluster, the stellar wind ejecta are completely re- 
tained and the surrounding medium remains unperturbed. The cur- 
rent mass of uj Cen is unlikely to generate sufficient dynamical fric- 
tion to modify its orbit to its present small size (Majewski et al. 
2000). Dynamical modeling, however, points out that the long and 
complex star formation history of wCen is inconsistent with the 
cluster originating on its present orbit: with a period of only 120 
Myr (Dinescu et al. 1999), the frequent disc crossings would have 
swept out all the intracluster gas very soon, leading to a mono- 



Figure 3. Comparison of observed (thin dashed histograms) and predicted 
(thick solid lines) metallicity distribution functions of uj Cen stars. Theoreti- 
cal predictions refer to the evolutive picture where uj Cen is the remnant of a 
larger system evolved in isolation and then accreted and partially disrupted 
by the Milky Way. The theoretical distributions are convolved with a Gaus- 
sian of dispersion <j = 0.2 dex in order to (generously) take observational 
errors into account. 

metallicity system. Thus, the progenitor of uj Cen must have been a 
massive enough system to allow dynamical friction to drag it to the 
inner Galactic regions (Bekki & Freeman 2003). In light of these 
considerations, in the following we discuss the results for a nucle- 
ated dwarf hosting uj Cen with initial mass ./#bar = 10 9 Mq . After a 
3 Gyr evolution, the parent system has lost most of its gas through 
galactic winds and ended up with a stellar mass .-# s tar S ~ 10 8 Mq, 
which is consistent with that of uj Cen's parent galaxy according to 
the computations of Bekki & Freeman (2003). 

4.2. 1 Stellar metallicity distribution and age-metallicity relation 

The observed MDF of long-lived stars in a galaxy is an important 
record of its past evolution. In fact, the relative numbers of stars 
which formed at any metallicity testify the interplay of fundamental 
processes such as star formation, infall of gas from the surround- 
ings and metal-enriched gas outflows at any time. Reproducing the 
currently observed MDF of uj Cen's stars significantly restrains the 
free parameter space of our dwarf galaxy model. We find that a fast 
early collapse coupled with an intense (per unit mass) star forma- 
tion activity, (ip) ~ 0. 1 Mq yr _1 during the first 1 Gyr evolution, 
gives rise to a distribution peaked at [Fe/H] ~ —1.6, as observed 
(Fig. 3, left panel), thanks to the strong galactic wind which effi- 
ciently removes the metals from the proto-oj Cen. In our model, the 
thermal energy of the gas exceeds its binding energy nearly 200 
million years after the onset of the star formation, owing to multi- 
ple SN explosions: the gas is then swept away by a strong galactic 
wind and the star formation slowly fades. We impose that the wind 
is differential, namely that the SN ejecta leave the galaxy more eas- 
ily than the stellar wind ejecta (see also Recchi et al. 2001). In 
particular, by assuming a wind efficiency Wi ~ 5 Gyr -1 for the 
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Figure 4. Predicted AMR for u> Cen (thick solid line) compared to (i) the 
AMR inferred from a sample of ~250 SGB stars with [Fe/H] errors less 
than 0.2 dex and age errors less than 2 Gyr, for two different choices of 
distance modulus and reddening (stars; Hilker et al. 2004); (ii) the age- 
metallicity diagram for the MSTO sample of Stanford et al. (2006; dots, 
only stars with V < 18 are considered). Notice, however, that the most 
metal-rich stars at [Fe/H] > —1.0 might have an accreted origin and ages 
comparable to that of the main cluster population (see text for references). 



SN ejecta, no stars with [Fe/H] > —0.4 dex are formed, in agree- 
ment with the observations (Fig. 3, left panel). Such high outflow 
rates - more than ten times the SFR - are often required to repro- 
duce the high-quality data of nearby dwarf spheroidals (Lanfranchi 
& Matteucci 2003, 2004). In Fig. 3 we compare the predictions of 
this model (thick solid lines) with the same observed distributions 
of Figs. 1 and 2. Both the steep rise at low metallicities and the ex- 
tended metallicity tail are well reproduced; in particular, the agree- 
ment with the up-to-date distribution of Sollima et al. (2005a) is 
strikingly good. This is encouraging, but the model has to be tested 
against many more observational constraints in order to prove its 
validity. 

In a recent study by Hilker et al. (2004; see also Kayser et al. 
2006), newly derived spectroscopic abundances of iron for ~400 
lu Cen members have been used in combination with the location 
of the stars in the CMD to infer the age-metallicity relation (AMR) 
of the system. The suggested age spread is about 3 Gyr and there 
is some indication that the AMR could level off above [Fe/H] ~ 
— 1.0 dex. This is clearly seen in Fig. 4, where the stars with the 
error bars represent the AMR for the subsample of stars with most 
reliable age and metallicity determinations, for two choices of the 
reddening and distance modulus values (see Hilker et al. 2004; their 
table 1). Evidence for the most metal-rich stellar populations of 
lu Cen being younger by 2-4 Gyr than the most metal-poor one has 
been found also by Stanford et al. (2006) from extensive Monte 
Carlo simulations on their MSTO data (Fig. 4, dots)§ . Although a 
large dispersion is present in the data, our model clearly predicts the 
correct run of metallicity with age (Fig. 4, thick solid line). Notice 
that in our model the SFR goes to zero after a 3 Gyr evolution, when 
[Fe/H] ~ —0.4 dex, because of the strong galactic outflows which 
efficiently remove any residual gas from the galaxy (see previous 
paragraph). 



3 Note that such an age difference seems inconsistent with the results ob- 
tained by Ferraro et al. (2004) and Sollima et al. (2005b) from the magni- 
tude level, shape and extension of the turnoff region. 
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Figure 5. Predicted (solid lines) versus observed (filled circles; see text 
for references) abundance ratios of several chemical species to iron as a 
function of [Fe/H]. The a elements oxygen (panel a), magnesium (panel 
b), silicon (panel c) and calcium (panel d) are shown, as well as the iron- 
peak element copper (panel e). Conservative errors are shown in the right 
hand corner of each panel. The crosses in panel e mark the time elapsed 
since the beginning of star formation (0.1, 0.3, 1, 2, 2.6 Gyr); they are not 
superimposed on the track in all panels in order to avoid confusion and 
overcrowded plots. 



4.2.2 Abundance ratios 

In Fig. 5 we display our predictions for several abundance ratios as 
a function of [Fe/H] (thick solid lines). The data displayed in Fig. 5 
(circles) have been collected from the literature. Data from high- 
resolution optical spectra (Francois, Spite & Spite 1988; Brown 
& Wallerstein 1993; Norris & Da Costa 1995; Smith et al. 1995, 
2000; Cunha et al. 2002; Pancino et al. 2002; Vanture et al. 2002) 
are shown as small circles. Data from low- and medium-resolution 
infrared spectra (Origlia et al. 2003) are shown as big circles. 
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We have checked that all literature abundance ratios are in rea- 
sonable agreement with each other (we consider two studies in rea- 
sonable agreement when the abundances of the stars in common do 
not differ, on average, by more than 0. 1 dex, which is the typical 
uncertainty of the measurements). As far as [Fe/H] is concerned, 
most papers agree with each other with the exceptions of Francois 
et al. (1988), Pancino et al. (2002), Vanture et al. (2002) and Origlia 
et al. (2003), all having [Fe/H] ~ 0.2 dex lower, which is a marginal 
(~ 2 <j) discrepancy. We therefore compared atmospheric parame- 
ters (Trff, log g and vt), atomic data (log gf) and equivalent width 
(EW) measurements among the above studies, to find possible ex- 
planations for the discrepancies. In the case of Francois et al., the 
turbolent velocities are significantly lower (~ 0.5 km s _1 ) and 
the adopted solar composition is significantly higher, log e(Fe)Q = 
7.67. Vanture et al. have instead a much lower v t (by 0.4 km s _1 ) 
than average, and log gf slightly lower (~ 0.06 dex). In the case 
of Pancino et al., all parameters appear in good agreement with the 
other literature sources, but for the only star in common (ROA 371) 
there is an average difference in EW of ~ 6.5 mA. All these factors 
together suggest that we should revise the Francois et al., Pancino 
et al. and Vanture et al. [Fe/H] values upwards by 0.2 dex. Origlia et 
al. have only stars in common with Pancino et al., with differences 
in T c ff (~ 150 K),log g (~ 0.3 dex) and v t (~ 0.3 kmT 1 ), but very 
similar abundances. Even if we do not have the tools to fully un- 
derstand the cause of the discrepancies (the method, resolution and 
spectral range employed by Origlia et al. are very different from the 
rest of the papers), we choose to revise the [Fe/H] values upwards 
by 0.2 dex, in conformity with what done for Francois et al. (1988), 
Pancino et al. (2002) and Vanture et al. (2002). 

Concerning the a-elements, there is a good agreement among 
various studies, except for a few notable exceptions. Oxygen is 
a very difficult element to measure, since only one line, [Ol] at 
6300 A, is used by most authors, that lies in a region plagued by 
atmospheric O2 absorption and is blended with a Ni I and a Sc II 
line. However, some authors also include the weak 6363 A line in 
their analysis, while only part of the authors perform a full spec- 
tral synthesis on the region. In spite of this, the abundance deter- 
minations are in reasonable agreement with each other, with the 
marginal exception of Vanture et al. (2002), that we choose to leave 
as it is, since the overall scatter in [O/Fe] is significantly larger than 
for other elements (see Fig. 5). For Mg I, Smith et al. (2000) have 
an abundance ~ 0.2 dex higher than average, which is partly ex- 
plained by the lower log gf (~ 0.1 dex) of the employed lines. We 
therefore lower the Smith et al. [Mg/Fe] abundances by 0.1 dex, 
bringing them in reasonable agreement with other determinations. 
Silicon abundances are a tricky business, since there is a general 
disagreement among authors (only Francois et al. 1988 and Smith 
et al. 2000 appear to be on the same scale), but no clear cause for 
the discrepancies emerges from our analysis, so we choose to apply 
no correction for this particular element. The Ca I determination by 
Norris & Da Costa (1995) is ~ 0.2 dex higher than all other studies 
considered here, that is fully taken into account by the ~ 0.2 dex 
lower log gf values adopted by those authors, so we revise Norris 
& Da Costa's [Ca/Fe] downwards by ~ 0.2 dex. 

Finally, only two papers have studied Cu in uCen so far, 
namely Pancino et al. (2002) and Cunha et al. (2002), having only 
one star in common, ROA 371. The linelist for the synthesis of the 
Cu I 5782 A line used is different in the two papers. The adopted 
atmospheric parameters are very similar, but the solar copper abun- 
dances differ by 0.09 dex. All this together makes an ~ 0.2 dex 
difference between the two studies of ROA 371 ([Cu/Fe] = —0.33 
in Pancino et al. and [Cu/Fe] = —0.52 in Cunha et al.) as can be 
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Figure 6. Theoretical Type II (solid line) and Type la (dashed line) SN rates 
(number per century) for the u> Cen progenitor. The SNIa rate is multiplied 
by a factor of ten to make it clearly visible. 

seen from fig. 6 of Cunha et al. (2002). However, if we adopt the 
Grevesse & Sauval (1998) solar composition for Cunha et al., as in 
Fig. 5, the discrepancy rises to 0.3 dex ([Cu/Fe] = —0.61 in Cunha 
et al.). We choose to use the same solar abundance (Grevesse & 
Sauval 1998), to put into evidence the real discrepancy between the 
two studies, which reflects the actual uncertainty in the derivation 
of Cu abundances. 

Theoretical values are normalized to the solar ones of 
Grevesse & Sauval (1998). Elements such as carbon and nitrogen 
are not considered, as their original abundances could be easily al- 
tered in the atmospheres of the sampled giant stars T The neutron 
capture elements Y, Ba, La and Eu deserve special attention and 
will be treated in detail elsewhere (Romano et al., in preparation). 

It can be seen that the model reproduces satisfactorily the 
decreasing trend of various [a/Fe] ratios with time (metallicity), 
with the possible exception of Mg, which is underestimated for 
[Fe/H] > —1.0. Appropriate adjustment of the adopted nucleosyn- 
thesis prescriptions and IMF slope would make the model predic- 
tions fit the data even better. Indeed, Francois et al. (2004) find that 
the Nomoto et al. (1997) yields of Mg we are using need signifi- 
cant corrections to best fit the abundance data of very metal-poor 
Galactic halo stars, and the stellar mass function is known to flat- 
ten in the very-low stellar mass domain (e.g. Kroupa et al. 1993; 
Chabrier 2003). However, what really matters here is to provide 
an overall, coherent interpretative framework for the whole body 
of data for ui Cen, rather than to attempt to precisely fit a specific 
observational constraint. 

The high ratios of [O/Fe], [Mg/Fe], [Si/Fe] and [Ca/Fe] for 
[Fe/H] < —1.2 clearly point to SNell as the major drivers of the 
chemical enrichment. Yet, the indication of a decrease of these 
abundance ratios in the intermediate and metal-rich subpopulations 
(Pancino et al. 2002; Origlia et al. 2003) are the unmistakable sign 
of significant Type la SN pollution, occuring at later times but be- 
fore star formation stops. 

As far as copper is concerned, it is worth stressing that the flat 
behaviour of the [Cu/Fe] data indicates no evolution in the copper- 



^ Actually, a large scatter is present also in the abundances of oxygen. 
Although commonly explained as the signature of processing in the CNO 
cycle (e.g. Norris & Da Costa 1995), it is worth stressing that it could be at 
least partly due to the wealth of observational problems discussed above. 
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to-iron yields over much of the chemical evolution within uj Cen 
(see Fig. 5, panel e, circles, and Cunha et al. 2002, their fig. 6). 
Yet, from Fig. 5, panel e, it is seen that our chemical evolution 
model (thick solid curve) overpredicts the copper-to-iron ratio over 
the whole metallicity range. The predicted behaviour of [Cu/Fe] 
results from the assumption that the major astrophysical site for 
the synthesis of Cu are SNela (Matteucci et al. 1993), and these 
stellar factories happen to contribute to the chemical enrichment 
of uj Cen already at the lowest metallicities. This is apparent from 
Fig. 6, where we show the Type II (solid line) and Type la (dashed 
line) SN rates predicted by our model for u> Cen. The SNIa rate has 
been multiplied by a factor of ten to make it clearly visible. SNell 
appear soon (a few million years) after their massive progenitors 
are born, thus their rate closely follows the SFR. SNela, instead, 
come from intermediate- to low-mass stellar progenitors in binary 
systems (Matteucci & Greggio 1986), which explode on varying 
time-scales: they start to contribute significantly to the chemical 
enrichment of the proto-tj Cen ~400 Myr after the beginning of 
star formation, when the ISM has attained a metallicity of [Fe/H] ~ 
— 1.6 (Fig. 6, upper x axis), while the maximum enrichment takes 

place only several Myr later, at t ~ 1.1 Gyr, when [Fe/H] ~ — 1 . 2 1 1 . 
Later on, at t ~ 2.6 Gyr ([Fe/H] ~ —0.6), the contribution from 
SNela nearly equals that from SNell, but only a few stars form from 
this SNIa-enriched material before the star formation stops. It is 
worth emphasizing here that the metal-enriched outflow lengthens 
the enrichment time-scale, i.e. the time that it takes for the ISM to 
reach a given metallicity: in the absence of such an enriched wind, 
it would take only 0.7 Gyr for the ISM to reach a metallicity of 
[Fe/H] ~ -1.2 and 1.8 Gyr to reach [Fe/H] ~ -0.6. 

In their 1993 paper, in order to fit the solar neighbourhood 
data, Matteucci and coworkers adopted Cu yields from SNela a 
factor of 100 higher than predicted by current SNIa models. When 
adopting the lower Cu yields from SNela predicted by Iwamoto 
et al. (1999), we find a flat behaviour of [Cu/Fe] versus [Fe/H] in 
uj Cen (Fig. 5, panel e, thin solid curve). While this is consistent 
with the data of Cunha et al. (2002), it does not account for the rise 
for [Fe/H] > - 1.0 pointed out by Pancino et al. (2002). Clearly, the 
issue of copper production in stars need to be further investigated. 
Apart from the flatness, an even more striking feature is the overall 
Cu deficiency in conjunction with a strong .s-process enhancement. 
The same chemical signature appears to characterize the Sagittar- 
ius dwarf spheroidal galaxy, thus lending support to the idea that 
uj Cen is the remaining nucleus of an accreted dwarf Galactic satel- 
lite (McWilliam & Smecker-Hane 2005). 



4.3 Helium enrichment in uj Cen 

In this section, we explore the issue of the helium enrichment in 
uj Cen within the picture where this cluster is the compact survivor 
of a larger satellite system ingested by the Milky Way many Gyr 
ago. 

In Fig. 7 we compare the theoretical AY versus [Fe/H] re- 
lation (thick solid line) to the data (box and filled circles). When 
using standard stellar yields (see Section 3.3), after a 3 Gyr evo- 
lution we obtain a negligible increase of the He abundance in the 
ISM (AY = 0.01), thus underestimating by an order of magnitude 



II Notice that a bimodal distribution of the delay times for the explosion, 
with 'prompt' and 'tardy' events, has recently proven to best match present 
SNIa data (Mannucci, Delia Valle & Panagia 2006). 
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Figure 7. Relative He abundance as a function of metallicity. The theoret- 
ical relation (thick line) is compared to the empirical one obtained from 
RR Lyra; star data (filled circles; see Sollima et al. 2006a, and references 
therein, for details about the derivation of the relative He abundances in 
these variable stars via the mass-luminosity parameter A). The box repre- 
sents the level of He enhancement required in order to explain the bMS data 
(Norris 2004; Piotto et al. 2005). 

the level of He enhancement required to explain the bMS and hor- 
izontal branch (HB) data (AY = 0.12-0.14; Norris 2004; Piotto 
et al. 2005). Yet, our model predictions are fully consistent with 
the tiny change in helium abundance (AY = 0.035 ± 0.066) over 
a wide metallicity range (—2.2 ^ [Fe/H] ^ —1.1) implied by the 
observations of 74 RR Lyrae variables by Sollima et al. (2006a). 
The fact that there is a population of metal-intermediate RR Lyras 
stars in u> Cen with luminosity and pulsational properties that are 
incompatible with a significant helium overabundance adds new 
complexity to the peculiar chemical features of uj Cen, as it en- 
tails the coexistence of two populations with similar metallicities 
but very different helium abundances in the cluster (Sollima et al. 
2006a). 

The stars on the blue side of the MS show no spread in metal- 
licity (Piotto et al. 2005), thus suggesting that they formed from a 
well-homogenized medium. On the other hand, in order to elevate 
Y from its primordial value 0.248 to ~0.40 one must assume that 
the material from which the bMS stars formed was made up al- 
most entirely of pure ejecta from previous stellar generations (Nor- 
ris 2004). In fact, as long as the fresh ejecta of dying stars are di- 
luted and mixed up with pre-existing gas in the framework of a ho- 
mogeneous chemical evolution model, there is no way of attaining 
a significant helium enrichment of the ISM by the time the metal- 
licity increases from [Fe/H] ~ —2.0 to ~ —1.2 dex (see Fig. 7, 
solid line). This result is independent of the adopted stellar yields, 
as can be immediately understood from an inspection of Figs. 8 and 
9. There, we show the mass fraction of He in the ejecta of LIMSs 
and massive stars as a function of stellar mass, calculated as 



- 1 ejecta — ^ mi 1 , V"/ 

m - m rem 

where Ymi is the initial helium abundance of the star, pHe is the stel- 
lar yield, namely the fractional mass of the star of initial mass m 
which is restored to the ISM in the form of newly produced He, 
and mrem is the mass of the remnant. The quantities displayed in 
Fig. 8 have been computed using the tables of van den Hoek & 
Groenewegen (1997) for LIMSs and Nomoto et al. (1997; stars) 
and Woosley & Weaver (1995; filled circles) for massive stars. In 
the case of van den Hoek & Groenewegen's yields, the bigger the 
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Figure 8. Mass fraction of He in the ejecta of LIMSs as well as massive 
stars as a function of stellar mass. V e jecta values were calculated using the 
tables by van den Hoek & Groenewegen (1997) for LIMSs and Woosley 
& Weaver (1995; filled circles) and Nomoto et al. (1997; stars) for massive 
stars. In the case of van den Hoek & Groenewegen's yields, the bigger the 
size of the symbol, the higher the initial metallicity of the star (Z = 0.001, 
0.004, 0.008). For comparison, the Y value suggested for the bMS stars of 
lu Cen is also displayed as a dashed line. The adopted stellar lifetimes are 
reported (for a few objects) on the upper x axis. 
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Figure 9. Mass fraction of He in the ejecta of rotating stars as a function 
of stellar mass. Y e jecta values were calculated using the tables by Meynet 
& Maeder (2002), Hirschi (2006) and Meynet et al. (2006) for different 
initial metallicities and/or rotational velocities of the stars. Upside-down 
triangles: Z = 10~ 8 ; stars: Z = 10~ 5 ; circles: Z = 0.004. Empty symbols: 
"ini = 800 km s -1 ; filled symbols: t>i n j = 300 km s —1 , except for the 9 Mq 
and 40 Mq stars at Z = 10~ 8 , for which v M = 500 km s" 1 and 700 km 
s _1 , respectively. For comparison, the Y value suggested for the bMS stars 
of uj Cen is also displayed as a dashed line. The adopted stellar lifetimes are 
reported (for a few objects) on the upper x axis. 



size of the symbol, the higher the initial metallicity of the star (Z = 
0.001, 0.004, 0.008). The quantities displayed in Fig. 9 have been 
computed using the tables of Meynet & Maeder (2002), Hirschi 
(2006) and Meynet et al. (2006) for rotating stars, for different ini- 
tial metallicities and/or rotational velocities of the stars (see fig- 
ure caption). We refer to those authors for details about the stellar 
model assumptions. In the case of the standard yields (Fig. 8), l^jecta 
is always below the quantity needed to explain the bMS data. Ro- 
tating stars, instead, do actually (slightly) exceed Y s j xla = 0.38 in 
the high-mass domain, but below 8 Mq they always contribute a 



much lower He amount (Fig. 9). Since, according to the Salpeter 
IMF we are using, in each stellar generation about 85 per cent of 
the mass falls in the 0. 1-8 Mq mass range, it is easy to guess that 
after the ~ 1 Gyr evolution needed to achieve [Fe/H] ~ — 1 .2 in the 
ISM, the average helium abundance of the medium out of which 
the bMS stars form will be definitely lower than required. Indeed, 
when including the yields of the Geneva group in our homogeneous 
chemical evolution model, we always find Y < 0.30 (AY < 0.05). 
We will come back to the issue of the helium enrichment in uj Cen 
in a future paper (Romano et al., in preparation). 



5 DISCUSSION AND CONCLUSIONS 

In this work we study the formation and evolution of the anomalous 
globular cluster uj Cen. The constraints established by the wealth 
of very good-quality abundance data for its member stars collected 
over the last decade significantly restrain the evolutive picture. We 
show that in the framework of the closed-box self-enrichment sce- 
nario (still an often used approximation) the metallicity distribution 
function of uj Cen's stars can not be reproduced. On the other hand, 
the main chemical properties of uj Cen are nicely reproduced if it is 
the compact remnant of a dwarf spheroidal galaxy evolved in iso- 
lation and then accreted - and partly disrupted - by the Milky Way. 
This evolutive picture has already allowed different authors to ex- 
plain the present-day peculiar location near the Galactic centre and 
surface brightness profile of uj Cen (see, e.g., Bekki & Freeman 
2003; Tsuchiya et al. 2004; Ideta & Makino 2004); our analysis 
gives further strength to the accreted dwarf picture by proving that 
the ingested satellite would also have similar chemical properties to 
uj Cen. Indeed, by assuming a relatively long-lasting star formation 
activity (though with most of the stars forming within 1 Gyr), stan- 
dard IMF and standard stellar yields, our models satisfactorily re- 
produce several observed abundance ratios as a function of [Fe/H]. 

Regarding the helium abundances of uj Cen, a great deal of 
work has been recently devoted to the enigma of the anomalous 
helium-to-metal enrichment of bMS and hot HB stars (Norris 2004; 
Lee et al. 2005; Piotto et al. 2005; Bekki & Norris 2006; Maeder & 
Meynet 2006). Stars which populate the 'blue hook' region could 
have suffered an unusually large mass loss on the RGB and, con- 
sequently, experienced the helium core flash while descending the 
white dwarf cooling curve (Moehler et al. 2002). Therefore, it is 
not clear which fraction of the He enhancement truly reflects the 
pristine abundance. Stars on the bMS, instead, likely give more 
trustworthy indication: they point to AY = 0.12-0.14 (AY/AZ ^ 
70). Recent observations of RR Lyrae stars suggest that a popula- 
tion with a normal (i.e. nearly primordial) helium content inhabits 
the cluster as well, thus complicating the overall picture (Sollima 
et al. 2006a). Our homogeneous chemical evolution model, when 
adopting a 3 Gyr long star formation history, standard yields and 
Salpeter's IMF, predicts that the helium abundance in the ISM al- 
most does not change during the cluster's progenitor evolution, in 
agreement with the RR Lyras star data, but in sharp disagreement 
with the abundances inferred for the bMS. We compare existing he- 
lium yields for low-metallicity rotating stars to the outputs of stan- 
dard stellar evolution which does not take rotation into account, and 
find that the winds of massive stars of intermediate rotation veloc- 
ities do indeed eject significant amounts of He (see also Maeder 

6 Meynet 2006). However, once weighted with a normal IMF, the 
helium yield of a stellar generation turns out to be lower than that 
required to explain the bMS anomalies (see Section 4.3). It is worth 
noticing at this point that even if the escape of SN ejecta through 
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galactic winds helps to keep the global metal content of the galaxy 
low, thus favouring high AY/AZ values, it also leads to longer 
enrichment time-scales (see Section 4.2.2), thus allowing low- and 
intermediate-mass stars to spread their relatively helium-poor gas 
through the ISM. Ad hoc scenarios in which the winds of massive 
stars remain confined and pollute only the central regions where 
the He-rich MS stars are actually more concentrated may represent 
a way out of the problem. Helium diffusion has recently been pro- 
posed as a viable mechanism to significantly increase the helium 
content of protostellar clouds (Chuzhoy 2006). 

Assuming an IMF strongly biased towards massive stars could 
also solve the problem, but current observational evidence seems to 
favour a Salpeter-like IMF in dwarf spheroidals as well as Galac- 
tic globulars (Wyse 2005, and references therein). Another option 
would be that the bMS was formed apart as a star cluster and then 
merged with the rMS in the central region of to Cen's host galaxy, 
a possibility first pointed out by Bekki & Norris (2006). We note 
that if the star cluster originated mostly from matter ejected by 
slow winds from rotating massive stars with [Fe/H] ~ —1.2, then 
it would have both the high helium abundance and uniform metal- 
licity inferred from bMS observations. Measuring the helium abun- 
dances of stars more metal-rich than [Fe/H] ~ — 1.0 would help us 
to discriminate between a scenario where the helium abundance in 
the ISM monotonically increases with time and a scenario where 
the formation of a helium-rich population is a fortuitous event oc- 
curred on a short time-scale under very special conditions. We will 
examine these scenarios in detail in a forthcoming paper, together 
with the most updated (even ad hoc) nucleosynthesis prescriptions. 
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